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Abstract 

The spin Green's function of the antiferromagnetic Heisenberg model on a trian- 
gular lattice is calculated using Mori's projection operator technique. At T = the 
spin excitation spectrum is shown to have gaps at the wave vectors of the classical 
Neel ordering. This points to the absence of the antiferromagnetic long-range order 
in the ground state. The calculated spin correlation on the neighboring sites of the 
same sublattice is in good agreement with the value derived from exact diagonal- 
ization. The temperature dependencies of the spin correlations and the gaps are 
calculated. 
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Strong electron correlations are expected to play an essential role in the recently discovered 
electron doped cobaltate Na x Co0 2 *?/H 2 [1,2,3,4]. This system has layered structure which 
consists of two-dimensional C0O2 planes separated by Na + and H 2 layers and resembles the 
structure of cuprate high-T c superconductors. However, in contrast to square Cu0 2 planes 
the Co ions form a triangular lattice. It is supposed that for moderate Na concentrations 
x the low-energy physics of Co0 2 planes can be described by the t-J model [5]. For small 
x the major part of Co ions are in the Co 4+ state with the spin S = 1/2 and the t-J 
model is reduced to the two-dimensional 5 = 1/2 Heisenberg model. The knowledge of the 
spin excitation spectrum of this model is necessary for further investigation of spectral and 
magnetic properties of the Co0 2 planes. Various numerical and analytical methods were 
used for elucidating the ground state properties of the Heisenberg model on a triangular 
lattice [6]. In Refs. [6,7] it is argued that the ground state of this system has the long-range 
antiferromagnetic order. The opposite conclusion was made in Refs. [8,9] analyzing the data 
of the exact diagonalization of small clusters. 
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In this work we use the method [10] based on Mori's projection operator technique [11] in 
which the spin Green's function is represented by a continued fraction. The elements of the 
fraction are calculated in the recursive procedure which is similar to Lanczos' orthogonal- 
ization. The decoupling is used in calculating the terms of the fraction. By analogy with 
Ref. [12] the decoupling is corrected by introducing vertex parameters the values of which 
are determined from the constraint of zero site magnetization following from the rotational 
symmetry of the Hamiltonian and the nearest-neighbor spin correlation derived from exact 
diagonalization. The used approach retains the rotation symmetry of spin components and 
does not contain any preset magnetic ordering. With this procedure the energies of mag- 
netic excitations and spin-spin correlations are calculated. The calculated spin correlation on 
neighboring sites of the same sublattice is in good agreement with the value available from 
exact diagonalization [8,9]. The spectrum of magnetic excitations is found to have finite gaps 
at wave vectors of the classical Neel ordering Q = (^£, , (^,0^ for zero temperature. 
This indicates that the ground state of the Heisenberg model on a triangular lattice has no 
long-range magnetic ordering. 

The Hamiltonian of this model reads 

H=lY,J™(< sZ m + 4 1 S m 1 ), (1) 



where s z n and s£ are the components of the spin-| operators, n and m label sites of the 
triangular lattice, a = ±1. We take into account nearest neighbor interactions only, J nm = 
J 12a ^n,m+a where J is the exchange constant and the six vectors a connect nearest neighbor 
sites. 

Magnetic properties of the model are inferred from the spin retarded Green's function 

D(kt) = - l 9(t)([st(t),s z _ k ]), (2) 



where s£ = N~ l l 2 ^ n exp(— ikn)s^, iV is the number of sites, s k (t) = exp(iHt)s k exp(—iHt) 
and (...) = Sp[exp(— H/T) . . .]/Sp[exp(— H/T)} with the temperature T. As mentioned, 
we apply Mori's projection operator technique [10,11]. To obtain the spin Green's function 
we use the continued fraction representation for Kubo's relaxation function ((<Skl s -k))w = 
/^dtexp(iu;OKI^k))i,M 

KKk))« = — s ' k V > ( 3 ) 

U) — E y~ 

oj - E x 



and the relation 

D(u;t)=u;((st\s z _ k %-(si,s\), (4) 
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where (A, B) = i J °° dt{[A(t), B}). 

The elements of the fraction and V* are determined from the recursive procedure 



[A n , H\ = E n A n + A n+1 + K_iA»-i, E n = (L4„, H], 4) 4)~\ 

(5) 

K-i = (A n , 4) (Ai-i) 4-0" 1 , ^-i = 0, A) = 4> n = 0, 1, 2, . . . 
The operators Aj constructed in this procedure form an orthogonal set, (A^Aj) oc 5^. 
Using procedure (5) we get 



k- 



E = («*, 5i k )(4, s^)" 1 = ([4, s* k ]>(4, si,)- 1 =0, A 1 = is 
F = 6J|d|(l -7 k )(4^- k ) _1 , #i = (i^.-w'JCwk.-i^k)" 1 = 0, 

where 7 k = ^cos(k x ) + | cos(4f) cos(^|^) in the orthogonal system of coordinates, C 1 = 
( s n ls n+a) * s the s P m correlation on neighboring sites. At this point we interrupt the continued 
fraction and calculate (s k , si k ). In the accepted approximation V\ oc (A 2 , A\) = 0. From this 
equation we find [10] 

(\i 2 st, -<^J> = 36^(1 - 7k ) 2 (4, ^k)" 1 . (6) 

The quantity i 2 s k in the left-hand side of this equation contains terms of the type sfs^s" 1 
which, following Refs. [12,13], are approximated by [(aiP nm + 0:2(1 — Pnm)){Sa l Sm)(^ ~ 
8nm) + |^nm]sf . Here ai t 2 are the vertex corrections, P nm = 1 if the sites n and m belong to 
the same magnetic sublattice and zero in the opposite case. In the decoupling the constraint 
of zero site magnetization 

«) = \ ~ (Sn 1 ^) = (7) 

was taken into account. This constraint is fulfilled both in the paramagnetic and in the 
ordered states; in the latter case due to the averaging over all possible directions of the 
magnetization. Using this approximation for i 2 s k from Eq. (6) we find (s k , si k ) and from 
Eqs. (3) and (4) we get 

D(M = «Ki^WI, (8 ) 



W — CO: 



k 



where 



* = - U ( A + \ + *) , A = -W^&V W 
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where C' 2 = |+2(s+ 1 s n j: a „) with the vector a" connecting nearest neighbor sites of a magnetic 
sublattice and C 2 = 2C\ + (s^s^a)- The quantity is the frequency of spin excitations 
and the parameter A describes the spin gap at the wave vectors Q of the classical Neel 
ordering. As seen from Eq. (9), the frequency of spin excitations tends to zero when k — > 
(the Goldstone mode) and has local minima at k = Q. 

To finish the calculations we have to find the parameters cxi, a 2 , C±, C 2 and C' 2 ' in Eq. (8) 
and (9). For this purpose we use Eq. (7) and the relations connecting the spin correlations 
with Green's function (8) 

= ^ X> (k (n - m))i^ coth (£) . (10) 

However, these relations - the three equations (10) for C±, C 2 , C 2 and the constraint (7) - 
give only four equations for the above five parameters. At zero temperature this deficiency 
can be made up by using results of exact diagonalization. More specifically, we use the value 
Ci = —0.1215 derived for an infinite crystal from results of exact diagonalization [8]. If the 
assumption is made that for T = the system has the long-range antiferromagnetic order 
then the frequencies of spin excitations at Q have to vanish. In this case one more unknown 
parameter - the condensation part C - appears in the calculations [13]. This part is the 
contribution of points Q to the spin correlations (10), 

) = C cos (Q (n - m)) + ?^ £ cos (k (n - m))^— ^. (11) 

Also the above set of equations is supplemented by the equation A = with A given by 
Eq. (9). The sublattice magnetization m is directly connected with the condensation part: 

m = \j3C/2. Solving the mentioned set of equations for an infinite crystal we found that 
C = —0.05. This result contradicts to the assumption of the long-range order for which C 
has to be positive. We arrived at the conclusion that either there is no long-range order 
in the ground state of the 2D Heisenberg model on the triangular lattice or the ordering 
parameter is too small to be determined by the present approximate method. Notice also 
that the value of (4 ls n+a") = 0-077 obtained in our calculations is in good agreement with 
the value 0.0803 ± 0.004 derived from exact diagonalization [8] (to our knowledge C\ and 
( s n ls n+a") are the only correlations with the values well-defined from exact diagonalization 
for the triangular lattice). One of the manifestations of the short-range order is a finite 
frequency ujq of spin excitations at the wave vectors Q. 

To consider the temperature variation of the spin correlations and spectrum we decreased the 
number of parameters, setting a± = a 2 = a. In this approximation we have four equations 
for four parameters. The range of temperatures < T < 1.3J and a 1000 x 1000 lattice with 
periodic boundary conditions were used in these calculations. The approximation ol\ = «2 
is cruder than the previous one where ot\ and 0:2 could differ. Therefore the gap uq in the 
spectrum of spin excitations at the wave vectors Q, the existence of which is an indication 
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of the short-range order, is somewhat larger in the former approximation than in the latter. 
The temperature dependencies of this gap, the spin correlation C\ and the dispersion of the 
spin excitations in the approximation a% = a 2 are shown in Fig. 1-3. The gap Uq grows 
with temperature which corresponds to the decrease of the correlation length. 

In summary, Mori's projection operator technique was used for investigating the excita- 
tion spectrum of the two-dimensional Heisenberg model on a triangular lattice. The zero- 
temperature spin-excitation spectrum has gaps at the wave vectors (^-, ^0, f^f,0j of the 
classical Neel ordering which is an indication of the lack of the long-range antiferromagnetic 
ordering in the ground state of the system. Temperature dependencies of the spin correlations 
and the gaps were calculated. 
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Figure captions 

Fig. 1. The dependence of the spin correlation C\ = (s„ 1 s~^ a ) on temperature. 

Fig. 2. The dependence of the energy of the spin excitations at the wave vector of the classical 
Neel ordering Q on temperature. 

Fig. 3. The dispersion of the spin excitations at zero temperature. 
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